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Abstract 

We present a method to implement 3-dimensional J\f=l SUSY Yang-Mills theory (a theory with 
two real supercharges containing gauge fields and an adjoint Majorana fermion) on the lattice, 
including a way to implement the Chern-Simons term present in this theory. At nonzero Chern- 
Simons number our implementation suffers from a sign problem which will make the numerical 
effort grow exponentially with volume. We also show that the theory with vanishing Chern-Simons 
number is anomalous; its partition function identically vanishes. 
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I. INTRODUCTION 



Lattice studies of supersymmetric field theories have long been an elusive goal. The lattice 
regularization generally breaks the supersymmetry and the IR limit of the lattice theory 
typically contains SUSY breaking fermion and scalar mass terms, which in 4 dimensions 
arise at all orders in perturbation theory. Much effort has recently been directed toward 
finding ways to avoid these terms by creative choices of lattice actions and/or generalized 
supersymmetry algebras [H, 0] . 

One way around this problem is to consider minimally supersymmetric Yang-Mills theory, 
which contains only gauge fields and fermions. The theory is then automatically (acciden- 
tally) supersymmetric provided one can correctly implement the fermions. Recent advances 
in fermion implementations [!, 0] have made it possible to achieve this program in 4 dimen- 
sions [EJ. 

In this paper we instead consider 3-dimensional minimally supersymmetric Yang-Mills 
theory (N=l supersymmetry, with two real supercharges). This theory is free of scalars 
and so correct implementation of the fermions again yields the right supersymmetric IR 
limit "accidentally." However, the implementation of the fermions is quite intricate, since 
one must impose a Majorana condition, and the implementation is further complicated by 
phases arising both from the fermionic determinant and from a Chern-Simons (CS) term, 
which is possible (and we will show, required) in this theory. 

The goal of this paper is to give a recipe for studying N—l SUSY in three dimensions 
with a Chern-Simons term on the lattice. This theory is believed to display very interesting 
nonperturbative properties that make it a prime target for simulation. In particular, Witten 
has conjectured [6| that the theory either preserves or spontaneously breaks supersymmetry, 
depending on the value of the Chern-Simons term. A lattice study of the theory would then 
constitute the first test (to our knowledge) of a nonperturbative supersymmetry breaking 
mechanism. 

The paper will be organized as follows. In section [Til we will review the continuum action 
of the theory, the necessity of a Chern-Simons term, and the anomaly condition which fixes 
the Chern-Simons term to take certain half-integer values. We will also present a proof, 
apparently unrecognized before, that the theory with vanishing Chern-Simons term has a 
vanishing partition function and is therefore not well defined. In section ITTT1 we will describe 
the discretization process, showing that the magnitude of the fermion determinant can be 
included using a rooted 3-D (2 component) Wilson-Dirac fermion with SW improvement and 
mass counterterm. We then show how to extract the Chern-Simons phase and the phase 
of the properly regularized rooted determinant. We leave the simulation itself for a future 
work; the goal here is to show that such simulations can be done and to provide the required 
tools. 



II. 3D N=l SYM IN THE CONTINUUM 



The field content of the theory consists of a gauge field and a 2-component, Majorana 
fermion in the adjoint representation (gaugino). In 4-dimensional notation, the gaugino 
is a 2-component Weyl fermion which has been further reduced from 2 complex to 2 real 
components by the application of a Majorana condition (possible in 3 dimensions). The 
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Euclidean action is 



s = j 2 j<Px {^f; v f; v + lum^J ■ (i) 

With these conventions the eigenvalues of the Dirac operator are all pure imaginary. Because 
the fermion is in a real representation of the gauge group, the fermionic operator possesses a 
doubled spectrum; if ip\ is an eigenvector of p with eigenvalue iX, then eip\ (where e = —IU2) 
is also an eigenvector with the same eigenvalue; 

o^D^x = iXipx 

ea^D^l = -ie\ip* x 
aMe^x) = iX(eiP* x ) ■ (2) 

The Majorana condition consists of taking only one of these degenerate sets of eigenvalues 
to define the fermion contribution to the path integral. That is, in the path integral replace 



detp = Y[Xi -> ^Jdetp = l[ , X l (3) 

i i 

where n/ is defined by taking only one eigenvalue from each degenerate pair. 

We can add a Chern-Simons term to this action in 3D provided we include an appropriate 
mass term for the fermions so that SUSY is retained: 

1 m r 7 g 2 k 
= -2itikN cs + -— %p%p , m= —— . 

2 g z Jx 47T 

Here k is the level of the CS theory or CS coupling. It is straightforward to check that the 
action is indeed invariant under the SUSY transformations 



i 
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and 5ip a = --e^ p F* u a p a, (5) 



with a the Grassman valued Majorana spinor parameterizing the SUSY transformation. 
Our convention for Majorana spinors is s = s T e. 

It has long been known 0, @] that - for gauge group SU(iV c ) - k must equal N c /2 modulo 
an integer to avoid a gauge anomaly so that, in particular, the theory with odd N c is ill 
defined for vanishing CS term. In Sec. Ill Al we review this argument and present a similarly 
motivated argument for the theory with even N c that implies that these theories are also ill 
defined for vanishing k. 



A. Anomalies in 3D SU(N C ) SYM 

In jsf Witten gave the first example of a theory with a rooted determinant that is sick 
with a global gauge anomaly, namely 4D SU(2) gauge theory with an odd number of left- 
handed fermion doublets. The problem with this theory is that it is impossible to define the 
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fermionic determinant over the space of gauge connections (gauge fields modulo all gauge 
transformations) such that it is both continuous and single-valued. 

We will review the 3D analog, formulated by Redlich [3] for gauge group SU(iV c ), which 
is slightly simpler and is of immediate importance for the current discussion. The general 
issue is that there are phase ambiguities in performing Grassman integrations. This is a 
problem in writing a path integral unless the phase ambiguity can be reduced to a single 
gauge-field independent phase, which factors out from the partition function and cancels 
in determining any correlation function. Therefore we pick some gauge field Aq, (chosen so 
that the Dirac operator has no zero eigenvalues) and call its contribution to the path integral 
real and positive. Then we determine the sign for any other configuration A by insisting on 
continuity of det p along a path from Aq to A. This is possible for any gauge field because 
our group manifold is path connected. 

The procedure is sketched in Fig. [TJ we watch the low lying eigenvalues of the spectrum 
as the path parameter t is varied from to 1 and count the number of eigenvalue pairs that 
change sign from positive to negative. This gives the relative sign of the two determinants 
in the path integral. 




FIG. 1: A path in configuration space and its associated eigenvalue flow 

This prescription is unique unless the number of eigenvalue zero-crossings depends on the 
path. This can happen because the space of connections is multiply connected; if two paths 
from Aq to A form a noncontractible loop, there is no guarantee that they lead to the same 
sign choice for detlp(A). 

Call the space of 3D gauge connections X. It contains noncontractible loops if the third 
homotopy group n 3 of the gauge group is nontrivial. (It is the third homotopy group because 
we are in 3 dimensions.) This is the case for all SU(iV c ), for which 7r 3 = Z. Consider a 
nontrivial loop from Aq back to Aq. This path is effectively a 4D gauge field configuration, 
where the space is S\X 3D space and the 4D gauge fields are A M = {A^.A^ = 0). The 
path is noncontractible if the instanton number of this gauge field configuration is nonzero. 
The four dimensional Weyl determinant has a number of zero modes determined by the 
Atiyah-Singer theorem [9J; for the fundamental representation this is 1 and for the adjoint 
representation this is 2N C . The 4D zeros correspond to zero crossings, and therefore sign 
flips, of the 3D fermionic determinant. If the number of sign flips in traversing the loop 
is odd, then the definition of the determinant cannot be both continuous, nontrivial, and 
single valued. 
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FIG. 2: A closed loop in configuration space with a non-trivial winding number 



For our case this is relevant because we want the square root of the Weyl determinant; 
the 2N C zero crossings become N c zero crossings when we choose one from each pair of 
eigenvalues, and this leads to a sign flip if N c is odd. 

There is an additional sign if the theory is defined with a nonzero Chern-Simons term. 
The CS term picks up a factor of 2nv in traversing a path of instanton number u, so the 
path integral picks up an overall factor of 

(-1)^ exp{i27Tuk) . (6) 

This implies that, in order to avoid a gauge anomaly, k=N c /2 modulo an integer. Only 
certain (half-integer) values of the Chern-Simons term are allowed, and in particular, the 
theory with N c odd and vanishing CS term is ill defined. 

We will now show that the supersymmetric theory with N c even and Chern-Simons 
coefficient k = (and therefore fermion mass of zero) is also anomalous, a point which 
to our knowledge has not been noticed before. Consider a configuration A(x) and its parity 
dual A'(x) = —A(—x). We claim that these give canceling contributions to the partition 
function if the Chern-Simons term is absent and the fermion mass is zero. They clearly have 
the same bosonic action, so we must show only that their rooted fermion determinants are 
equal and opposite. We define the sign of the rooted determinant for configuration A(x) to 
be positive and connect them with a path from A(x) to the trivial vacuum and from the 
trivial vacuum to A'(x) via the parity dual of this path (see again Fig. [[]). 

The Dirac operator for the trivial vacuum configuration has — 1) pairs of zero eigen- 
values (we implicitly work on a torus with standard boundary conditions). There are then 
— 1) pairs of eigenvalues that cross zero at the vacuum configuration. Furthermore, if 
n + pairs cross zero somewhere on the path between A(x) and the vacuum, than the number 
of pairs n_ that cross zero between the vacuum and A'(x) will be the same. The total 
number of pairs which change sign in going from A(x) to A'{x) is therefore (iV c 2 — 1) + 2n + . 
This is odd, so the fermion rooted determinant flips sign and the configurations give cancel- 
ing contributions to the partition function, which vanishes identically. An example of the 
eigenvalue flow for gauge group SU(2) is shown in Fig. [3] for the case m=0. 

This problem is avoided at nonzero Chern-Simons term because A(x) and A'{x) have op- 
posite Chern-Simons number and so enter the partition function with opposite phase, rather 
than canceling. Similarly, at nonzero fermion mass the eigenvalues of the Dirac operator are 
complex and introduce nonzero phases which are opposite between the configuration and its 
parity dual. 
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A(x) Vacuum A '« 




FIG. 3: Flow of almost zero mode eigenvalues of the Dirac operator as the configuration is varied 
between parity conjugate gauge fields on a path through the vacuum. 




Re A [lattice units] 

FIG. 4: Eigenvalues for 20 Wilson fermion configurations in a 8 3 box at g 2 a = 0.5, illustrating the 
"return loops" in the complex plane. 

B. Regularization dependence of determinant 

It is necessary to clarify what is meant by the "continuum" fermionic determinant. The 
issue is that continuum theories are always denned as limits of discrete (or otherwise regular- 
ized) theories and in a discrete theory with volume regularization there are a finite number 
of eigenvalues for the Dirac equation. In traversing a closed loop with nonzero instanton 
number, we just saw that a nonzero number of eigenvalues cross from negative to positive 
value. Since the final configuration is the same as the starting one, it has the same spectrum. 
Since there are a finite number of negative eigenvalues, there must be some compensating 
flow of positive eigenvalues to negative, somewhere in the complex plane. In other words, 
eigenvalues must "return" somewhere in the complex plane. We illustrate this for Wilson 
fermions in Fig. HI The figure superimposes the eigenvalue spectra, in the complex plane, of 
20 quenched gauge field configurations in an 8 3 box with lattice spacing g 2 a = 0.5. Each dot 
is a pair of eigenvalues (the pairing of the spectrum discussed in the last section occurs for 
both the Wilson and overlap lattice implementations of the Dirac operator). The spectrum 
of eigenvalues for p parallels the imaginary axis near zero but bends out into the complex 
plane for large eigenvalues A ~ 1/a and forms a loop, so eigenvalues moving from negative 
to positive values "push" eigenvalues around the loop to reappear at negative values. These 
extra vanishing-imaginary-part eigenvalues are the lattice fermion doublers which have been 
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pushed out into the complex plane by the Wilson term; for the Wilson action in 3 dimensions 
there are actually 3 extra places where eigenvalues cross zero imaginary part, corresponding 
to (n, 0,0), (ir, 7r,0), and (7r,7T7r) type doublers. 

The problem is that this "return loop" will contribute to the partition function even 
in the continuum limit. Because it represents deeply UV physics, its contribution to the 
partition function must be representable in an effective IR description in terms of local 
effective operators. The 3D theory has only one marginal operator, which is the Chern- 
Simons term. Therefore the return loop can (and will) induce a Chern-Simons term, but 
will not otherwise change the infrared description (in the small a limit). It is easy to see 
that the size of this Chern-Simons term is fixed by the rate of spectral flow near zero; the 
number of eigenvalues which circle around the loop must be the number needed to refill the 
negative eigenvalues when the spectrum flows upwards. The sign of this extra Chern-Simons 
term depends on whether the return loop is at positive or negative imaginary part, which is 
a regularization detail. The desired continuum theory is the one in which this regularization 
effect has been removed by a Chern-Simons number counterterm. 

III. THE DISCRETIZATION 

This section describes the discretization of N=l SYM with a Chern-Simons term in 3D. 
The theory contains nontrivial phases; that is, different gauge field configurations contribute 
to the partition function with different complex phase as well as different magnitude. The 
plan is to treat this using the Edinburgh method; one studies the theory on the lattice by 



building a Markov chain sample weighted by the magnitude of the action | exp(— S)y det Jp\ 
and then includes the phase as part of the observable. Phase cancellation reduces the 
statistical power in a volume dependent way. Therefore our implementation is on the same 
footing as finite chemical potential simulations in QCD; they work in principle, but whether 
they work in practice depends on how severe the phase cancellation problem turns out to 
be. 

The implementation consists of two parts; the real part of the bosonic action and mag- 
nitude of the fermionic determinant, and the phase. The real bosonic action is completely 
standard. We describe the magnitude of the determinant first, then the Chern-Simons part 
of the phase, then the phase in the determinant. 

A. Fermion Implementation 

Though there may be some advantages to using an overlap fermion implementation to do 
the simulation, we believe that the numerical simplicity of the Wilson implementation makes 
it a much more sensible choice. The usual problems with Wilson fermions are less severe 
than in 4 dimensions; chiral symmetry is a non-issue because there is no such thing in 3 
dimensions, and the additive renormalization of the mass is well behaved because the theory 
is super-renormalizable. One can easily work at a fine lattice spacing where the problem 
of exceptional configurations is well under control (note that we are only interested in the 
theory at finite fermion mass, since as we just argued the massless theory is anomalous). 

The fermionic action reads 





X 
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■ A „t; ( r - cr^U^(x)^ x + (r + <r ll )Uj l (x—n)ip as - lt 



2a 



(7) 



We can improve the convergence of the spectrum to the continuum limit from 0(a) to 0(a 2 ) 
via the 3D analog of the Sheikholeslami-Wohlert term [l(| (SW) term 



S sw = a 3 ^7rVU°>> a A ( p ^{x) - PL(x))ij; : 



X, /i, V 
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(8) 



Here P^ v (x) is the average of the 4 plaquettes in the \w plane as shown in Fig. [5j This 
improvement is probably necessary to implement simulations at reasonable lattice spacings. 



Uu 




^Tr (ul(x)UUx+fi)U^x+v)U„(x) + ... 



FIG. 5: Clover field strength for the Sheikholeslami-Wohlert improvement term. 



Since the 3D theory is superrenormalizable we need only the tree level determination of 
the SW coefficient, c sw =l, to remove 0(a) corrections other than mass renormalizations. 
Much of the complication of the improvement in 4D is thus avoided in the 3D version. 

The improvement to the spectral properties is quite dramatic, though it parallels closely 



the improvement in 4D so we refer the interested reader to, for example, [llj for analysis 
of 4D Dirac operator spectra with improved and unimproved actions in a situation which is 
fairly analogous to the cases we consider in 3D. Fig. [6] shows an example of the improvement 
to the physical branch of the spectrum for 20 configurations on a 14 3 lattice with g 2 a = 
0.5. Improvement "squeezes" the eigenvalues toward the solid line-a parabola incorporating 
the "bending" in the complex plane present in the tree level Dirac operator due to the k 
dependence of the Wilson term. The spectra shown are shifted appropriately so that both 
represent a fermion determinant for the case of zero physical mass. 

The disadvantage of Wilson fermions is that we must tune the fermion mass to remove an 
additive correction. To consistently improve the theory to eliminate all 0(a) errors we need 
to do this at the two loop level (see [12J for a discussion). Here we present the calculation 
of the one-loop mass counterterm; in practice both one and two loop counterterms can 
be determined quite easily numerically by analyzing the low-lying eigenvalues in the Dirac 
operator spectrum at a few different lattice spacings. 



The SW term modifies only the i/jAi/j vertex. Following closely the 4D treatment of 13 
we determine the new three point vertex to be (with ip(p), 4>(p') an d A(k=p' —p)) 



be 



la^ cos ■ 



(p+p') 



+ 



[It, a 



a I COS — 

n 2 



(9) 
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Re \ [lattice units] Re A [lattice units] 

FIG. 6: Physical branch of the spectrum of the Wilson-Dirac operator; 20 superimposed configu- 
rations (left) with and (right) without improvement (14 3 lattices with g 2 a = 0.5) 



The one- loop bare mass counterterm required to tune the theory to its SUSY limit, for 
Wilson coefficient r = 1, is 

5m = ^^(-2.3260) = -^(-2.3260) , (10) 

4% llXp 

with Ca the first Casimir of the adjoint representation of the gauge group - Ca=N c for 
SU(AQ. 

To eliminate all 0(a) errors it is also necessary to compute a 1-loop fermionic wave- 
function renormalization (which can be interpreted as an 0(a) multiplicative renormalization 
of the mass) and a 1-loop renormalization of the field strength (or equivalently of g 2 , or 
equivalently of the lattice spacing). The contribution of bosonic loops to the gauge coupling 
renormalization was computed in [13], and the fermionic contribution is not difficult. In the 
same notation as that previous work, 

Aatt(x) = E Tr C 1 - P ^ X )) + Z ^ i D w + D sw + f^m + ^] V , (11) 

Z g9 a i^j V z i> Z i>) 

where it was determined that 

vbosonic = l _^_( ^L(2N 2 C - 3) + N c f^Zi _ lX) W ith f = 0.152859325 , (12) 
9 in 1 9N r v c ' 1 12 9 ) ) S ' K ' 



we have 



?P- = 1- ^-^e m with e m = 2.75066732(8) (13) 
4tx 

and Zj otal = Z 9 bosonic - ^^e g with e g = 0.204254488(7) . (14) 
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To treat Nf fundamental representation fermions rather than one adjoint with a Majorana 
condition, replace N c — > Cf in the expression for Z m jZ$ and change N c — > Nf in the 
expression for Z g (the 1/2 trace normalization group theory factor cancels a factor of 2 from 
not imposing the Majorana condition). 

As already mentioned, 8m is easily measured nonperturbatively by - for example - looking 
at the lowest lying eigenvalues of the Dirac operator for a range of (3 values and fitting to find 
what value of 5m is needed to get the zero imaginary part eigenvalues to have vanishing real 
part as well. To tune out the 0(a) error analytically would require a 2-loop determination 
of the fermion self energy, which is much more difficult. 

The numerical implementation of the magnitude of the determinant is achieved by includ- 
ing (det P^T') 1 / 4 in the path integral, which can be accomplished by conventional pseudofer- 
mion techniques. There should be no issues with locality (we believe) because the operator's 
spectrum is doubled, so the fourth root is actually well defined (up to a phase). 



B. Evaluation of Chern-Simons number 

Here we detail how to determine the Chern-Simons number of a 3-dimensional lattice 
configuration. Ou r ap proach is borrowed from an investigation of the electroweak sphaleron 



rate by one of us [15|, but we summarize it here for completeness. 



The literal definition of Chern-Simons number Ncs i n terms of the integral of F A A — 
|A A A A A, given in Eq. (j3J), is too gauge dependent and lattice spacing sensitive to be 
of much use, so we use instead the following (gauge invariant) properties of Chern-Simons 
number, which can be taken as a definition of Nqs, modulo an integer: 

• iVcs for the vacuum is (modulo an integer); and 

• the N cs difference between two configurations is (modulo an integer) equal to the 
integral of e^ va ^ Tr F tlu F a f3/32TT 2 along a path through configuration space connecting 
those configurations. 

The second point requires some clarification. As previously discussed, given two 3D configu- 
rations A\(pc), A|(x), one can find a path connecting them through gauge field configuration 
space; that is, one can find v4*(x, r) with r £ [0, 1] an affine parameter and A l (x, 0) = ^(x), 
v4*(x, 1) = A|(x). We may think of the path as a 4-dimensional gauge field configuration, 
with t as the fourth coordinate, D = D T , and F 0i = i[D , Dj\. The Nqs difference between 
two configurations is the integral 

N CS (A 2 ) - = J d 3 x £ dre l]k ^ ^f ik ■ (15) 

We have not specified how A is to be chosen along the r direction, but this turns out not 
to matter; a different choice leads to a change in F 0i of Di5A , but the Di can be integrated 
by parts onto Fjj. and e^DiF^ vanishes by the Bianchi identity. The idea is then to define 
iVcs for a configuration (modulo an integer) by finding a path from that configuration to 
the vacuum and integrating FF along that path. 

The complication in using this approach to define iVcs on the lattice is that there is no 
lattice definition of the field strength F^ v which satisfies the Bianchi identity; therefore the 
procedure is ambiguous. Further, no specific prescription for choosing Aq generically leads 



to an integer iVcs around a closed loop. There are mathematically rigorous [16|, |17|, llSj and 
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numerically implementable [l^, 20] methods to find the integer value around closed loops, but 
these are not helpful here, since we really want the value on a path with distinct beginning 
and end points. The trick instead is to note that the problems with lattice implementations 
of FF arise when the fields are "coarse" (plaquettes far from identity, most of excitations at 
the lattice spacing scale). We can define iV cs uniquely by choosing any unique prescription 
for a path from a configuration to the vacuum. We can ensure that the result is as close as 
possible to the continuum meaning of Nq$ if our unique prescription is one which quickly 
smooths out the lattice-spacing scale fluctuations in the gauge field configuration. The early 
"smoothing out" part of the path then contributes a small UV "lattice artifact" contribution 
to iVcs and the remaining path gives a contribution which closely resembles the continuum 
value of iVcs for this configuration. 

A good choice is the "cooling path" or gradient decent under the energy, 

D rA*,r) = - 9 q%£^ ] , H[A} = /A~Tr4L4(y,r)]. (16) 



We use a Symanzik [21| or "rectangle- improved" definition of H and of the field strength 
appearing in FF, as described in THf, which also shows extensive tests of the approach. 



To summarize, the procedure is to determine the Chern-Simons number of a lattice 
configuration by integrating FF along the "cooling" or gradient-decent path through con- 
figurations to the vacuum. The procedure is unique and gives an answer which is continuous 
over the space of gauge field configurations except at a "sphaleron" separatrix, where it 
is discontinuous by (almost exactly) an integer. If Markov-chain configurations are tightly 
enough sampled, one can determine the integer part by continuity. 

On fine lattices this definition of Chern-Simons number should correctly reproduce the 
continuum notion up to corrections suppressed by two powers of the lattice spacing. 

A straightforward alternative method to implement integer Chern-Simons numbers would 
be to evaluate the phase in the determinant of a fundamental representation Wilson fermion 
with a negative mass of order the lattice spacing (so the origin of the complex plane is inside 
the leftmost "circle" in Fig. H]). This method would be numerically much less efficient, since 
the numerical effort in taking a determinant rises as the third power of the number of lattice 
points or a 9 , while the approach presented only grows worse as a 5 and can be reduced to a 3 
through the careful use of blocking [TEJ . 



C. Fermionic phase 

Since in practice only the magnitude of the rooted fermion determinant can be included 
dynamically in a lattice simulation, we must still describe a prescription for assigning to 
the Wilson-Dirac operator a phase for each field configuration in order to complete the 
prescription for the discretization of the theory. 

One approach to do this would be to explicitly evaluate the (complex) determinant of 
T> for each configuration in the Markov chain and then halve the phase, fixing the sign 
ambiguity by using a tightly sampled Markov chain and demanding continuity. Then one 
must subtract the Chern-Simons phase contribution from the "return loop" discussed in 
Subsection III Bl This approach is correct but numerically expensive. 

After correcting for the Chern-Simons term induced by the "return loop" this phase 
is dominated by the contributions of the low lying eigenvalues, so we can develop a more 
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efficient procedure by focusing on determining the phase arising from these eigenvalues. Since 
this point is key to our approach we should explain it in a little detail. Look again at Fig. [61 
In the continuum limit the eigenvalues will lie, not on an arc, but on a straight line with real 
part m and imaginary part set by the eigenvalue under the Tj) operator (for the free theory, by 
k). At large eigenvalue the weak-coupling approximation is valid (since the theory is super- 
renormalizable this is true by a power of the eigenvalue). The density of eigenvalues therefore 
scales with the free theory density of states, k 2 dk. An eigenvalue's phase difference from ±^ 
is tan _1 (m//c), so naively the phase arising from large eigenvalues could be large. What is 
important, though, is the phase difference configuration by configuration, and this becomes 
small, essentially because the large eigenvalues do not change very much as a function of 
the gauge field. To see this, note that the large eigenvalues represent short-range physics. 
The influence on the effective IR behavior can be expanded as Wilsonian renormalization 
of effective IR operators. The lowest order parity-odd operator is the Chern-Simons term; 
all others are higher by at least two powers of derivatives and therefore suppressed by at 
least g 4 /k 2 . Since we are allowing 0(a 2 ) errors, the only operator we need to incorporate 
correctly from the large k eigenvalues is therefore their contribution to the Chern-Simons 
term. Therefore our strategy will be to include low eigenvalues' phases explicitly and to 
determine the phase contribution of large eigenvalues in terms of their contribution to an 
effective Chern-Simons term. 

We can extract the smallest M eigenvalues using the Arnoldi method at much less nu- 
merical effort than is required to determine the full determinant. Each such eigenvalue \ 
of the Wilson-Dirac operator with renormalized mass m takes a value in the complex plane, 
(rj,Zj). In the continuum limit the real parts rj always equal m. At finite lattice spacing 
there will be 0(g A a 2 ) and 0(i 2 a 2 ) deviations in the real part. (There would be 0(g 2 a) devi- 
ations arising from the dimension 5 operator ipa^F^ip had we not included the SW term.) 
We "clean" the low lying eigenvalues by projecting them to the r = m axis, see Fig. [7J Each 
eigenvalue then contributes a phase 

(f>, = tan- ( & " m)2 + g) 
\ m J 

This amounts to projecting the physical branch of the Wilson-Dirac spectrum onto the axis 
of the continuum spectrum and then calculating the phase as sketched in Fig. [7J The idea is 
to incorporate the phase of all eigenvalues for which the angle fa lies in a range [— max , ma x]- 
In the a — > limit we must take max — > f ■ This requires more eigenvalues at finer lattice 
spacing; this can be made more efficient by using the shifted Arnoldi method. 

The Markov evolution of the gauge fields will move around the eigenvalues so that 
eigenvalues regularly move in and out of the "window" in which we explicitly include them. 
When an eigenvalue goes above or below max , the phase we determine will abruptly change 
by T^max- Therefore each configuration along the Markov chain must be reasonably close 
to the last, so that these phases can be determined by continuity. The change in phase 
between neighboring configurations C/i )2 contributed by all eigenvalues lying above max o.nd 
below — max is simply an effective Chern-Simons term, as discussed above. The size of the 
contribution can be determined by considering the amount of spectral flow due to a changing 
Chern-Simons number, and is well approximated by 2(| — (f) max )N c (Ncs(U2) — Nqs(Ui)) . We 
have confirmed this in quenched simulations, for instance by analyzing the max dependence 
of the procedure and seeing that this contribution ensures independence on this artificial 

12 



(17) 



FIG. 7: Sketch of the projection for the phase prescription. 



parameter. If we choose the vacuum with, for example, iVcs = to have zero phase, then 
this prescription uniquely determines a phase for all configurations. 



D. The sign (or phase) problem 

Now that the appropriate implementations of the fermionic field content and the Chern- 
Simons term have been detailed there remains no fundamental barrier to the simulation of 
the theory and so we turn our attention to an important technical issue of the simulation. 

The lattice simulation consists of replacing the path integral by a sum over a finite set of 
link field configurations that are distributed with a probability given by the Boltzmann factor 
for the theory, | exp(— £?[?/]) |. For complex action, the ensemble average for an observable 

will be 

M»«, (18) 

where = argexp(— 5" [[/]). Obviously this leads to cancellations between link configura- 
tions from the sampling, and thus to a reduction in statistics. That is, for a sample of N 
independent configurations, the error in the numerator scales as \/~N ; but the denominator 
will be smaller than N, so the error in the operator will not be 1/ \/jV. This problem could 
be eliminated by performing "phase quenching" on the theory, but this does uncontrolled 
damage to the theory which in our case we believe is severe. Therefore we must face this 
phase cancellation issue. 

To determine how bad the phase cancellations will be, we start by looking at the theory 
with very large fermion mass, so that the effect of integrating out the fermion is well approxi- 
mated by a shift to the Chern-Simons coefficient, k — > k — N c /2. In large volumes, we expect 
that N cs will be Gaussian distributed around zero. The degree of phase cancellation is deter- 
mined by how badly our determination of 1 in the denominator of Eq. ( 1181) is "suppressed" : 
all measurables must be scaled by the result of the partition function which we evaluate as 
(1) with the average replaced by a sum on configurations with phases. Fluctuations go as 

1 / \/N with iV the number of configurations (unless each configuration contributes better 



13 



statistics-this depends on the measurable). The point is that these fluctuations are to be 
compared with an average value which suffers phase cancellation. We estimate this by doing 
the (Gaussian) integral over N cs to find how big the partition function actually is: 



(1) 



dN r . 



expH27rfciV cs ) exp(-A 2 /2 X ) 



(19) 



with x the variance of N cs and k the effective Chern-Simons coefficient. The integral gives 

exp(-27r 2 k 2 X ) . (20) 
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FIG. 8: (color online) Variance of Ncs as a function of volume, for two lattice spacings, in quenched 
SU(2) gauge theory and using the definition of Nqq, presented in the text. 

Fig. [8]shows the dependence of \jV on the volume of the lattice, in quenched simulations. 
As the figure shows, going to a finer lattice induces a constant shift in x, and the volume 
dependence becomes weak for boxes larger than about 8/g 2 . At this value, on a "reasonably 
fine" lattice of g 2 a = 0.5 [(3 = 8], the Acs variance is about .175, leading to a sign problem 
induced loss of statistical power of order exp(— 2n 2 k 2 x .175) = 1/30 for k — 1. This volume 
is therefore at the limit of practicality. Statistical power falls exponentially if we try to 
make the volume any larger. However, a box of length 8/g 2 is very large. The theory has 
a mass gap and the large volume limit should be approached rapidly in a box a few times 
longer than the longest correlation length. For the SU(2) theory, the lightest glueball mass is 
mn++ ~ 1.66g 2 [22| and the inverse correlation length involved in the Debye mass is ~ 1.14g 2 
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These both suggest that the dominant physics is on quite short scales ~ l/g 2 , though 
the string tension suggests a longer correlation length ~ 3/g 2 Therefore this 

volume is probably sufficient to effectively achieve the continuum limit. 

The sign problem grows more severe at large k, as indicated by Eq. ( 1201) . Fortunately 
the mass scales with fc, so we can reduce the volume as V oc 1/k 3 in the large k limit. We 
must also tighten the lattice spacing a to keep ka fixed, and since (Nq S ) turns out to have a 
linear UV divergence (which causes the lattice spacing dependence observed in Fig. [8]), this 
means that x wm scale as y/a ~ k~ 2 . In the large k limit the severity of the sign problem 
therefore approaches a finite limit. 

Fig. [8] is based on quenched configurations. When we include the effects of dynamical 
fermions, we expect the situation at small k to improve for two reasons. The first is that 
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inclusion of the magnitude of the determinant in the Boltzmann factor suppresses configura- 
tions with large iVcs (an effect we have observed in preliminary quenched simulations). This 
is not surprising since we know that Dirac operators for configurations with half-integer iVcs 
(sphalerons) have zero eigenvalues. We thus expect this suppression to be more effective at 
smaller fermion masses. We will not attempt to quantify this statement further. 

The second reason is that the phase contributed by the determinant partially cancels 
against the phase from the Chern-Simons term. For large fermion mass, for example, we 
know from continuum methods that integrating out the fermion gives a shift in k or, equiv- 
alently, a contribution to the partition function exp(iirN c Ncs)- Furthermore, we have con- 
firmed for the case of N c =2 Wilson-Dirac fermions - with the prescription described above 
for assigning determinant phases - that this is approximately the phase of the rooted fer- 
mion determinant for the relatively small masses (k=l, 2, 3) of interest here. An example is 
shown in Fig. Ofor k—1. In this case the cancellation of the phase in the partition function 
is nearly exact, and the phase problem disappears. Finally we remark that, in the algorithm 
for determining Nqs, one can begin the integration of FF after a short amount of cooling, 
eliminating the contribution from the most UV modes. This eliminates a lattice artifact 
"noise" contribution to x but leaves that from interesting IR physics intact. 
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FIG. 9: Evolution of the phase of the rooted fermion determinant over a Markov chain of length 
200. The continuous line is 2ttN cs . The system transitions through a sphaleron and subsequently 
fluctuates around the topologically distinct vacuum at 2tt. 



IV. CONCLUSION 

Three dimensional minimally supersymmetric gauge theory can be implemented on the 
lattice. As we showed, this theory necessarily involves complex phases; the theory with 
vanishing Chern-Simons term and massless fermions is anomalous because the path integral 
is odd under parity transformations and the partition function therefore vanishes identically. 

Our implementation has concentrated on numerical efficiency Lattice spacing corrections 
should first occur at 0(a 2 ). Fermions need only be implemented using the clover-improved 
Wilson method, not the much more expensive overlap method. The Chern-Simons phase 
can be determined without reference to fermionic operators, and the phase in the rooted 
Dirac determinant is determined by finding only the low lying eigenvalues using the Arnoldi 
method, rather than by taking the full determinant. This efficiency is important because 
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the theory suffers from a sign problem which will make it difficult to take the large volume 
limit. We are guardedly optimistic that the sign problem will not be as severe as might be 
feared. First, the phase of the fermionic determinant and of the Chern-Simons term are of 
opposite sign and partially cancel. Second, when the Chern-Simons term is large, the theory 
is massive and the volume requirement should be reduced. 

Now that an implementable method has been presented, it would be very interesting 
to study 3D Af—1 SYM on the lattice. In particular it would greatly improve our insight 
into nonperturbative supersymmetry breaking if we could study and (presumably) verify 
Witten's conjectures regarding spontaneous SUSY breaking in this theory 
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